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ON ■ In both the large N c limit and the valence quark model, semileptonic de- 

ON 

' cays are dominated by resonant final states. Using Bjorken's sum rule in an 

a: 

O ■ "unquenched" version of the quark model, I demonstrate that in the heavy 

■ quark limit nonresonant final states should also be produced at a significant 



rate. By calculating the individual strengths of a large number of exclusive 



5h 

■ two-body nonresonant channels, I show that the total rate for such processes 

is highly fragmented. I also describe some very substantial duality-violating 
suppression factors which reduce the inclusive nonresonant rate to a few per- 
cent of the total semileptonic rate for the finite quark masses of B decay, and 
comment on the importance of nonresonant decays as testing grounds for very 
basic ideas on the structure, strength, and significance of the qq sea and on 
quark-hadron duality in QCD. 
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I. INTRODUCTION 



A. Background 

In both the large N c limit |lj and the valence quark model fHH, semileptonic heavy 
quark decays are saturated by resonant final states. In nature this idealization is broken by 
light quark pair creation which gives these infinitely narrow resonances widths and allows 
nonresonant final states to appear. While it is clear that qq pairs play an important role in 
this and many other phenomena, it is also clear that they remain poorly understood: 

1. Although from their widths it is easy to show that hadrons are full of qq pairs, meson 
and baryon spectroscopies are characterized by the valence degrees of freedom. In 
particular there is no evidence for excitations of the qq sea with respect to the valence 
quarks. 

2. Even if one were to assume that the qq sea is frozen out of spectroscopy, it is easy 
to show that, unless there is a conspiracy between them, meson loop graphs ought to 
destroy the successes of quark model spectroscopy. For example, the relative shift of 
the p and a% due to meson loops formed from their dominant decay modes (nir and 
pit, respectively) is hundreds of MeV H. 

3. Related to 2), but even more dramatic, is the relative shift of pairs of mesons such 
as the p and u. They are degenerate in the quark model and in the large N c limit, 
as observed in nature, but meson loop diagrams would lead one to expect them to 
develop a mass difference of hundreds of MeV || . 

4. Given that meson loop diagrams which generate the qq sea are very strong, it is sur- 
prising that the valence quarks seem to dominate low energy current matrix elements 
[[J. The extreme interest generated by the proton spin crisis may be attributed to the 
fact that it indicates that there are some current matrix elements where the valence 
quarks do not dominate. 
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While the questions qq pair creation raises are ubiquitous, heavy quark systems (e.g., 
Qd and udQ) are likely to be the most favorable systems in which to find answers: of all 
systems governed by strong QCD J?|, they are arguably the simplest. Indeed, many of their 
properties can be rigorously derived directly from QCD using heavy quark symmetry [^,0. 
Their simplicity also makes them more amenable to modelling than other hadrons. E.g., as 
"the hydrogen atoms of QCD" with the heavy quark defining an origin of coordinates, the 
simplest relativistic constituent quark models can treat Qd systems using the Dirac equation 
[TqH rather than a Bethe-Salpeter-type equation. Because the heavy quark is removed from 



consideration as mg — ► oo, such systems offer unique opportunities to study the "brown 
muck" one chunk at a time. 

In this paper I expand the quark model treatment of the "brown muck" in heavy meson 
semileptonic decay from a simple valence d or u antiquark confined to Q to include the 
leading effects of qq pair creation. See Fig. 1. Previous studies have examined low energy 
pion emission in the context of heavy quark chiral perturbation theory this work is 



the first of which I am aware to address the full array of nonresonant processes. While the 
results presented here will be specific to heavy quark systems, I will draw lessons from them 
and suggest experimental consequences of much wider interest. 

B. Nonresonant Final States in B and D Semileptonic Decay 

For D decays induced by the underlying c — > siv^ quark decay, D — ► Kiue and D — > 
K*lvi are clearly dominant, accounting for more than 90% of the inclusive D — > X s £v£ 
semileptonic rate fll2 |. For B decays induced by b — > civ?, B — > Dlvi and B — > D*£i>£ 



account for 64±6% of the inclusive B — > X c £v£ semileptonic rate [12|. Given the large energy 



release in b — > c versus c — > s, it is not surprising that in the former the two ground state 
resonances of the s[ l = | multiplet would account for less of the total semileptonic rate. 
Indeed, it has been argued |13[ by assuming duality between the ISGW2 valence quark model 
and QCD-corrected inclusive calculations jnj that a complete ISGW2-based calculation 
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would predict that 20 ± 6% of the B — > X c £z/^ rate should go to resonant excited states 
above the D and D*. (The quoted error corresponds to an estimate of the theoretical 
uncertainties in the QCD-corrected inclusive rate calculations.) ISGW2 as published is in 
contrast not exhaustive: it explicitly computes the rates to the L — 1 excited states with 
s'J e = \ + and s'1 l = | + and to the first radial excitations with s'J 1 = \ . These six 
lowest-lying excitations give 8 ± 1% of the QCD-corrected inclusive rate, implying that an 
additional 12 ± 6% of the rate of a complete the ISGW2 calculation should be in yet more 
highly excited states (both ordinary mesons and hybrids). The ISGW2 model and next-to- 
leading-order QCD can be compared in this way because both are valence-quark-plus-glue 
calculations: qq pairs are ignored in ISGW2 as 1/N C corrections and would enter the partonic 
level inclusive rates only at order a 2 s {via 1 /^-suppressed graphs). 



(a) 



(O 



AAA 



AA/ 1 



1 + r\A/v 



Fig. 1: Leading corrections to narrow resonance saturation of the rates for semileptonic 
heavy quark decay, with quark level diagrams on the left and their hadronic counterparts 
on the right. Internal hadronic lines are to be understood as summed over their full valence 
spectra, a) The valence — > valence graph with its leading vertex and external leg correc- 
tions, b) The valence — > valence graph followed by final state decay, and c) Decay from a 
nonvalence component of the initial state to a nonresonant two particle final state. 
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Experimentally, the extent of resonance dominance of B and D semileptonic decays 
remains unclear. In D decay there are explicit measurements giving D — > KkIvi rates 
which are 3 ± 1% of the semileptonic rate. Here the signal is excluded from being the 
K*, but it is not excluded that it could arise from the tails of broad resonances, so even this 
small fraction cannot be unambiguously identified as nonresonant. In B decay, 36 ± 6% of 
the semileptonic rate is not in the D or D* . Since the dynamical part p 2 dyn = p 2 — | of the 
slope of the Isgur- Wise function £ (it?) could be as large as about | fl5| , the loss of rate from 
these channels relative to p 2 dyn = in the approximation £ ~ 1 — p 2 {w — 1) could indeed be as 
large as 36%, i.e., the observed non- (£) + £)*) rate is not inconsistent with that expected from 
Bjorken's sum rule |jl6| , |17| in the heavy quark limit. Since, as just explained, one expects 
20 ± 6% of this non-(.D + D*) rate to be in excited resonant states, there would be room for 
16 ± 8% of B semileptonic decays to be nonresonant, i.e., for a large fraction of p 2 dyn to be 
due to nonresonant channels. There is another closely related indication that nonresonant 
channels might be significant: ISGW2 overpredicts the D and D* rates by amounts which 
are consistent with the observation that the Isgur- Wise function is falling about 25% more 
rapidly than expected from the opening of just excited resonance decay channels. While 
suggestive, both of these observations are also consistent with the ISGW2 model simply 



underpredicting decay rates to excited charm states |18|1 . An additional concern is that 
the suggestive loss of rate from the D and D* channels calculated from Bjorken's sum rule 
is actually an upper bound in the heavy quark limit: expected quadratic terms in £ will 
dampen this loss, suggesting that perhaps the experimental non-(D + D*) fraction is high. 

In summary, there is weak circumstantial evidence for nonresonant processes in heavy 
quark semileptonic decays. The most compelling case for the existence of such processes at 
some level is nevertheless the simple theoretical observation that in the real world we expect 
very strong nonvalence Qqqd components in a Qd state and so expect some inclusive rate 
to be lost from resonances and transferred to continua. In what follows I will make these 
qualitative expectations, previously outlined in Ref. |TB[, concrete. 
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C. The ISGW and ISGW2 Models 



In addition to providing a useful phenomenological guide to semileptonic decays, the 
ISGW model [0] was in many respects a stepping-stone to heavy quark symmetry, as it 
manifested this symmetry near zero recoil. ISGW2 [[| is an update of ISGW with many new 
features required by heavy quark symmetry: it includes constraints on the relations between 
form factors away from zero recoil and on the slopes of form factors near zero recoil [IT), 17], 



it relates the naive currents of the quark model to the full weak currents via the matching 
conditions of Heavy Quark Effective Theory (HQET) [JIT]], and it modifies the ad hoc ISGW 
prescription for connecting quark model form factors to physical form factors to be consis- 
tent with the constraints of heavy-quark-symmetry-breaking at order I/tuq. Several other 
improvements were also made, including the addition of heavy-quark-symmetry-breaking 
color magnetic interactions to the quark model's dynamics, the incorporation of relativistic 
corrections to the axial coupling constants (known to be important in the analogous coupling 
qa in neutron beta decay), and the use of more realistic form factor shapes, based on the 
measured pion form factor. For a more complete discussion of the foundations, strengths, 
and weaknesses of such models, see Refs. @,[|. 

In this paper I focus on correcting these models' neglect of the nonvalence components 
of hadrons, but note that many of the results to be described here are anticipated in the 
extensive discussion of Ref. [[| leading to the conclusion that nonresonant states could be 
ignored as a first approximation to the dynamics of semileptonic decays. 

D. Unquenching the Quark Model 

Some of the key puzzles associated with the nature and importance of qq pairs in low 
energy hadron structure were described above. These puzzles and potential solutions to 
them have been extensively discussed in a series of papers on "unquenching" the quark 
model [HH3]. In the following I briefly summarize these solutions, since they are the basis 
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for the study presented here. 



1. The Origin and Resiliency of Potential Models 

A central puzzle in hadron spectroscopy is the apparent absence of low energy degrees 
of freedom beyond those which can be attributed to the valence quarks (e.g., gluonic or 
sea quark excitations). Very closely related to this puzzle is the apparent unimportance of 
strong meson loop corrections. 

A simple resolution of this puzzle has been proposed f20 |. In the flux tube model ||2T 



the quark potential model arises from an adiabatic approximation to the gluonic and extra 
qq degrees of freedom embodied in the flux tube. This physics has an analog at short 
distances where perturbation theory applies. There Nf types of light qq pairs shift (in lowest 
order) the coefficient of the Coulombic potential from a^(Q 2 ) = „„. j^frnr; to ai Nt \Q 2 ) = 



33to(<2 2 1 Ag) 

12tt 



(33-2 jv f )[ n (Q 2 /A i .) ' ^he ne ^ e ^ ec ^ °f sucn pahs thus being to produce a new effective short 



distance QQ potential. Similarly, when pairs bubble up in the flux tube (i.e., when the flux 
tube breaks to create a Qq plus qQ system and then "heals" back to QQ), their net effect is 
to cause a shift AE^ f (r) in the ground state gluonic energy which in turn produces a new 
long-range effective QQ potential. 

It has indeed been shown [|J that the net long-distance effect of the bubbles is to create 
a new string tension b N (i.e., that the potential remains linear). Since this string tension is 
to be associated with the observed string tension, after renormalization pair creation has no 
effect on the long-distance structure of the quark model in the adiabatic approximation. Thus 
the net effect of mass shifts from pair creation is much smaller than one would naively expect 
from the magnitude of typical hadronic widths: such shifts can only arise from nonadiabatic 
effects [g2j. 

It should be emphasized that no simple truncation of the set of all meson loop graphs can 
reproduce such results: to recover the adiabatic approximation requires summing over large 
towers of Qq plus qQ intermediate states to saturate their duality with qq loop diagrams 
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which have strength at high energy. 



2. The Survival of the OZI Rule 

There is another puzzle of hadronic dynamics which is reminiscent of this one: the success 
of the OZI rule A generic OZI- violating amplitude Aozi can be shown to vanish like 
1/N C , and this is often quoted as a rationale for the OZI rule. However, there are several 



unsatisfactory features of this "solution" to the OZI mixing problem |24 ]. Consider u-(fi 
mixing as an example. This mixing receives a contribution from the virtual hadronic loop 
process u —>■ KK —>■ <fi, both steps of which are OZI-allowed, and each of which scales with 
N c like T 1 / 2 ~ iV" 1 / 2 . The large N c result that this OZI-violating amplitude behaves like 
N^ 1 is thus not peculiar to large N c : it just arises from "unitarity" in the sense that the 
real and imaginary parts of a generic hadronic loop diagram will have the same dependence 

on N c . The usual interpretation of the OZI rule in this case that "double hairpin 

graphs" are dramatically suppressed is untenable in the light of these OZI-allowed 

loop diagrams. They expose the deficiency of the large N c argument since Aozi ~ T is not 
a good representation of the OZI rule. (Continuing to use u-ip mixing as an example, we 
note that m w — is numerically comparable to a typical hadronic width, so the large N c 
result would predict an u-cj) mixing angle of order unity in contrast to the observed pattern 
of very weak mixing which implies that Aozi « T << m.) 

Unquenching the quark model thus endangers the naive quark model's agreement with 
the OZI rule. It has been shown || how this disaster is naturally averted in the flux tube 
model through a "miraculous" set of cancellations between mesonic loop diagrams consisting 
of apparently unrelated sets of mesons (e.g., the KK, KK* + K*K, and K*K* loops tend 
to strongly cancel against loops containing a K or K* plus one of the four strange mesons of 
the L = 1 meson nonets). Of course the "miracle" occurs for a good reason: the sum of all 



hadronic loops is dual to a closed qq loop created and destroyed by a Pq operator [25,26 



but in the closure approximation such an operator cannot create mixing in other than a 



S 



scalar channel. It can also be shown || that current matrix elements like S'y^s vanish in 
this same approximation. 

3. A Summary Comment on Modelling the Effects of qq Pairs 

The preceding discussion strongly suggests that models which have not addressed the 
effects of unquenching on spectroscopy and the OZI rule should be viewed very skeptically 
as models of the effects of the qq sea on hadron structure: large towers of mesonic loops 
are required to understand how quarkonium spectroscopy and the OZI rule survive once 
strong pair creation is turned on. In particular, while pion and kaon loops (which tend to 
break the closure approximation due to their exceptional masses) have a special role to play, 
they will not allow a satisfactory solution to these fundamental problems associated with 
unquenching the quark model and so cannot be expected to provide a reliable guide to the 
physics of qq pairs. 

Indeed, I hope the reader can appreciate just on the basis of this lightning review that 
there are great dangers in drawing conclusions about the strength, structure, or significance 
of qq pairs in hadrons from any model that has not dealt with these issues. 

II. UNQUENCHING HEAVY QUARK DECAY 
A. Background 

To unquench predictions of the quark model for semileptonic heavy quark decay, I will 
apply without alteration the model of Refs. [^-0 which solves the phenomenological prob- 
lems associated with unquenching the quark model. In particular, I assume that the qq 
pair is created by the action of a pair creation Hamiltonian density H pc (x) in the Qd flux 
tube. I further assume that the pair is created with a nonlocality (corresponding to a finite 
constituent quark radius) in the coordinate v. See Fig. 2. 
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Fig. 2: The coordinates for a) the Qd system, and b) the Qqqd system. In each diagram 
the cross x denotes the location of the center-of-mass; most results presented in this paper 
are in the heavy quark limit where the center-of-mass coincides with the position of Q. 

The coordinates used here in Qd are the standard center-of-mass and relative coordinates 

mQd 

r = r~ d -r Q (2) 

while in Qqqd the choice is 

^ = m Q?Q + m q (r q + fg) + m d r 3 ^ ^ 

m Qqqd 

r = r d --r Q (4) 

v = r q -r q (5) 

u> = -(r, + r ff ) ~ -(r g + r q -) -r Q , (6) 

where = m « + + ... + m^. (Note that the three- vector coordinate w should not be 

confused with the Lorentz invariant heavy quark scalar product w = v' ■ v called "double u" 
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(or, in many European countries, "double v") used in heavy quark form factors). Inverting, 
we have in Qd 

r Q = Rem ~ e d/Q d r ~ R cm (7) 

Td = Rem + tQ/Qd ? ~ Rem + T (8) 

where e a /p = m a /mp. In Qqqd we have 

?Q = R'cm - tqq/Qqqd ™ ~ 6 d/Qd ? ~ R'cm (9) 

fg = R' cm + e Q d /Qqq d w+^v~ R'^ + w + -v (10) 

1-1 

r q = R' cm + eQd/Qggd w--v^ Km + w--v (11) 

rd = R'cm - £qq/Qqqd W + € Q/Q g f ~ R' cm + f . (12) 

Note that most of the results of this paper are presented in the heavy quark limit where the 
approximations shown in these formulas will often be used. 
Thus for a CT state 

*Q- = 7^e iP --%d(f )x° SQSs (13) 
where X° SQ s d * s the s P m zero wavefunction, so that 

is the momentum space Qd wavefunction with 

<P Q M = - V / d're-^Qdir ) (15) 

and accordingly 

\PQd(Pcm)) = y/2m QS J d 3 p<p Q d(p )x° SQScl \Q( e Q/QdPcm - P, s Q )d{e d/Qd P cm + p, sj)> . (16) 



Note that in the limit P cm — > relevant here, the factor ^j2m Qd ~ ^2E Pq - is purely 
conventional as itlq — > oo. 

When the flux-tube-breaking pair creation Hamiltonian = / d 3 xH^(0, x) acts, 
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Hpcl-PQd(-Pcm)) — Vgg\PQqqd(Pcm)) (17) 

where, according to the flux tube model, 

1. since the qq pair is created in the flux tube, its center-of-mass is found in a wavefunction 
ipft(w,r) defined by the flux tube spatial profile, 

2. the internal wavefunction of the qq pair has J PC = ++ and, independent of w and f, 
is of the form ?p™ c (v) ■ X^ q where 4>^(v) has L = 1, Xj^ q nas S = 1, and ip m ■ x~ m = 
-^(^X- 1 - ip°x° + Y), and 

3. the amplitude to find the Qd subsystem inside the Qqqd system at relative separation 
r is identical to that in the ground state, namely i^ngif ), since the pair creation 
Hamiltonian density acts locally and instantaneously on the flux tube. 

In this formulation, r\ qq - defines the strength of the pair creation, and the normalized state 
\PQq q d(Pcm)) is determined by the wavefunction just described: 

*W = t^j^^M*, ? )<W ■ x^odf )x° SQSd (is) 

The component parts of this wavefunction and of Eq. are defined by previous studies. 
From ISGW2 [|, 

n3/2 

^Qd-^tji e 2 Qd (19) 

where /3qj = 0.41 GeV as ttlq ->ooas determined variationally from a Coulomb-plus-linear- 
plus-hyp erfine Schrodinger equation. The pair creation wavefunction ip pc is constrained in 
Ref. [|J by fitting decay data assuming the form 

J P c ~ ve- 3v2/Sr « = ve-^> 2 (20) 

to have a quark radius < r q < 0.4 fm. Given this constraint, I will take r q = 0.3 fm as 
a "canonical" value, corresponding to (3 pc ~ 0.58 GeV, but will consider deviations of ±0.1 
fm from this value as plausible. (This central value and range are guided by the difficulty of 
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inventing a mechanism which could lead to a constituent quark radius r q < 0.2 fm.) Ideally 
26| , ifjf t (w,r ) should have a probability profile which is a tube around the Qd axis with 



"caps" at Q and d. This structure is probably very significant for the decays of highly 
excited states, but since all our decays will emerge from the Qd ground state, I adopt a 
simpler and more heuristic model which simply takes 

„3/2 

= (21) 
with (3f t = fb, b being the string tension and / a coefficient with "canonical" value 2 and 



an uncertainty estimated to be ±1 based on the calculations of the properties of the ground 
state wavefunction of the flux tube presented in Appendix A of Ref. [p6fl . Some defects 
in this simplification and some subtleties associated with both the alternative of using a 
flux-tube shape and nonrelativistic kinematics are discussed in Appendix A of this paper. 

Since strong decay amplitudes are determined by matrix elements of H^(t, x) between 
the decaying particle and the continuum, r] qq can be determined empirically. In Appendix 
B, I extract from D* — > Dtt and K* — > Kn decays the coefficients r\ qq for q = u or d. For 
concreteness, I will assume following Refs. |2£],§J, that r] sS is identical, but that the t]qq for 
Q = c,b,t are all zero. Finally, for ease of exposition, I will treat explicitly the case of a 
single qq flavor in what follows, but take into account q — u, d, s in numerical results. 

Within these approximations, 

= s 3 (P - P m )M#)<l%(iO ■ x; q %^(p )x° SQSd (22) 

where (f>f t , 0™, and 0qj are the momentum space wavefunctions corresponding to i/jf t , ii™, 
and ipQd, respectively, so that 

\P Q q q d{Pcm)) = prn Q - qqS J d 3 u; J d^ J d 3 p ■ X7 q %<P Q( i(P )x° SQS3 

\Q{ e Q/Qqqd~Pcm ~ e Q/Qd^ ~ V-, s Q)^.{ e qj QqqdPcm + ~3 + ^ ■, Sq) 
<l{ e q/QqqdPcm + ^ ~ 7T, S q ) d(t d / Qqq ^P cm — e d / Q ^LO + f), Sj)) . (23) 

(Up to this point, I have retained the exact kinematics of the nonrelativistic limit for finite 
rriQ, but from now on I will generally simplify results by taking the heavy quark symmetry 
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— # — * 

limit nig — > oo with V cm = P C m/ m Q fixed.) 

While Eq. (|T|]) defines the action of H^J on the Qd sector, it does not of course provide 
us with the H^-perturbed ground state. This state is of the form 

\Pq{Pcw)) = 7= \PQd{Pcm)) + Cqq\pQqqd(Pcm)) (24) 

^/ 1+C2 qq 

where \pQq q d(P C m)) is a normalized Qqqd state of the same general form as Eq. (^) but 
with a wavefunction $Qg g( j to be specified below and determined via 

d 3 q 



C qq\PQqqd{Pcm)) — / 



aft 17 2m Qqqd 



\P cm ;ab(q))(P cm ;ab(q)\U^\P QS (P cm )) 



E a b(P C m, q ) — Ep Qd (P cm ) 

Here |P cm ;a6(g )) is the two meson eigenstate with Qq in internal state a, qd in internal 
state b, and with (Qq) a and (gcfJb having relative momentum <f and total momentum P cm ; 
E a b(Pcm, q ) and Ep Qd (P cm ) are the total energies of the respective Qqqd and Q<i states 
with fixed total center-of-mass momentum P cm . To obtain a rough expression for I-Pq^j), I 
make use of the approximate duality between each of the towers of states a and b and their 
corresponding free particle spectra in the internal relative momenta pqq and p q d, respectively. 
E.g., I use 

Z)((Q?)a| ~12 I d3 PQq(Q(-PQqi S Q)q(PQqi S q)\ ( 26 ) 

a SQSq J 

where for simplicity I have illustrated the duality equation in the Qq center of mass frame. 
While this replacement is imperfect for low a (corresponding to small PQq), since is 
quite pointlike, Eq. Q25D has most of its strength for relatively massive states a and b and 
the use of duality should be satisfactory for our purposes. 

With these approximations we can change variables from a, b, q to pQ q , p q g, q, and then 
convert both sides of Eq. (p5|) to the variables uj, fr, and p to identify 

r $ _ ^ Vqq^Qqqd / 97 N 
°qq^Qqqd ~ l^'J 
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where in the rest frame P cm = and with itlq — > oo, 

2 2 

= 2m g + ^— + + 5 (28) 

with 



m q Am q 



p 2 

5 = m Q + m d + — - m p QS . (29) 

In the duality approximation we have adopted, and with the wavefunction of the Qd system 
identical in Qqqd and Qd, 6 ~ 0. Moreover, using the wavefunctions fl2"0|) and (pID and the 
parameters given earlier 



< 4m 9 > 3/3^ 



< -i > 20/52 



(30) 



is small so that 



TV 2 



AE ~ 2m g + — (31) 



m 9 



and we may deduce that 

^6 S (P- P cn Mft^W P l c ^) ■ xZ l M(v )X° SQS , (32) 

i.e., that $Qg 9 j differs from &Qq q d simply by the replacement 4> pc — > (f> pc where 

s ^tWU ■ (33) 



The normalization factor n is given by 



n=[*«\ W*> J 2 ■ (34) 

which may be quite well approximated by the formula n = 1+21 1 ^ 12 ^ 3 with y = From 
Eqs. (p7|) and (|33|) follows the key relation that 



Note that </>^(vr) is softer than 0™(7r). We will often use a harmonic approximation 
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to (f)™. It is of course the softer shape of 0™ that makes n < 1, and thus (3 pc can be 
determined by requiring that Eqs. (|33| ) and ( p^) match near 7f = 0, i.e., that 

P pc = n 1 ' 5 P pc . (37) 

With realistic parameter values, 0.7, so the softening is not dramatic: with our 

canonical value for @ pc , [3 pc ~ 0.4 GeV. 

B. The Unquenched Isgur-Wise Function 

We are now in a position to calculate the unquenched quark model contribution to the 
Isgur-Wise function ||. By heavy quark symmetry, the form factors for a general Qi — > Q2 
transition can be calculated as matrix elements for the simpler Q —> Q transition with an 
arbitrary current QTQ. We therefore focus on the matrix elements of the scalar current QQ 
between Q-containing states: 

f °" w = ^ {PQ{ ^T )mPQ{ -^f )) (38) 



where w = v' ■ v ~ 1 + P c 2 m /2mg = 1 + and where, 

C H = ^(^-(^)IWI^-(-^)) (40) 

= / d 3 rr QS (r y- imd ^ T ^ Q - d (r ) (41) 

and 

— * — * 

<-QM i„,.\ — 1 1 td 7^imlP. .i'^ 



= — (W-F)iwiW~F)> (42) 



egf(^)er^) (43) 



with 
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& M (w) = / d 3 w^%(w)e- 2im « v ° m - a Tp ft (w) . (44) 



(In these equations the notation "QM" reminds us that we are calculating in the quark 
model so that the perturbative matching of these HQET matrix elements to field theory 
must be done at the quark model scale \lqm ~ 1 GeV.) 

We see from Eq. that does not depend on the poorly known qq wavefunction 

ip pc . This simplification arises because this qq wave function defines the relative position 
of the q and q, while £f t M is sensitive only to the qq system's wave function relative to Q. 
Defining for w — 1 << 1 

C QM ^l-p 2 wf (w-l) (45) 

C - 1 - PU W - 1) ( 46 ) 
~ 1 1) (47) 

and recalling the conventional definition 

£~l-p 2 (w-l) (48) 

we therefore have (displaying now explicitly the effects of summing over q = u, d, and s) 
simply 

T c 2 o 2 

n 2 - n 2 4- q J3LM (AQ\ 

Pwf ~ pQd + 1 7 v o (49) 

L ' J->q u qq 

where p 2 w ^ is the nonrelativistic wavefunction contribution to p 2 , i.e., it excludes the rela- 
tivistic | and the contribution Ap 2 t from matching to the low energy effective theory 0. 
Using Eqs. (f|l]) and (|43f) we then have the old result 

2 m 2 < r 2 > m 2 

^ 5 ^ (50) 



and the new correction from qq pairs 
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so that in the SU(3) limit where m u = rrid = m s = m 



2 2m - 



2 rn 2 d ^ c w-^f m[ d M 

- tJ- + APL • (53) 

This is one of our main new results. It shows that even if the c 2 = are large, the contribution 
of pairs to p 2 will be small in the adiabatic limit where they are highly localized in the flux 
tube (i.e., as 0j t — > oo). See Appendix A for a discussion of this result for a more general 
flux-tube shape. We will see next that to the extent that pairs contribute to the exclusive 
"elastic" slope p 2 , they will also contribute to the inclusive nonresonant semileptonic rate. 



C. A Duality Interpretation of Ap 2 ea via Bjorken's Sum Rule 

1. Motivation 

We have just seen that even if Pq x is full of qq pairs, they may not contribute to p 2 . We 
will now see that it is incorrect to associate the relative probabilities of Q\d and Qiqqd in 
Pq 1 with the resonant and nonresonant parts, respectively, of Qi —>■ Q2 semileptonic decay. 
As heavy quark symmetry requires, at w — 1 the Q\ — > Q2 transition creates only Pq 2 
and Vq 2 of the ground state s'J e — ~ multiplet independent of the structure of the "brown 
muck". I.e., in this limit the Q2d and Q2qqd components of the hadronic final state, no 
matter what their relative strengths, form perfectly into the resonant states Pq 2 and Vq 2 . 
For w — 1 small but nonzero, Bjorken's sum rule |TB,T7] tells us that the loss of rate from 



the "elastic" transitions Pq 1 — ► Pq 2 and Pq 1 — > Vq 2 relative to structureless hadrons with 
p 2 = 1/4 will be exactly compensated by the production of s'J 1 = \ + and | + states. In the 



valence quark model, this rate must appear in Q 2 d excited states. In Ref. |L7fl this valence 
quark model duality to the quark level semileptonic decay was explicitly demonstrated. Here 
I will show that the Qiqqd content of Pq x leads in general to the production of both resonant 
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and nonresonant final states, with the latter rates proportional to Ap^ ea . In particular, in 
the adiabatic limit there will be no nonresonant production. 

To the extent that Ap 2 sea is nonzero, nonresonant (Q2q)a(qd)b final states with s'J 1 = \ + 
and | + will be produced to compensate for the additional loss of rate from the elastic 
channels which it causes. It is natural to expect the compensation to occur in these channels. 
The softening of the elastic form factors which depletes the rate to Pq 2 and Vq 2 will have its 
analog in inelastic resonance excitation form factors, so these rates will also be diminished 
and cannot compensate for the additional loss of rate from Pq 2 and Vq 2 . The population 
of inelastic channels is thus the only avenue available for satisfying Bjorken's sum rule. Of 
course this is also intuitively appealing: the loss of rate to the elastic channels occurs because 
after the recoil from —P cm /2 to +P cm /2, the qq parts of the ground state wave functions 
of the initial and final states fail to overlap, and in so doing they must "by conservation 
of probability" find themselves in their excited states, namely as {Q2<l)a( ( ld)b continua. We 
will now make these heuristic observations precise. 



2. Production of Nonresonant States at Low w — 1 



The n th valence state perturbed by (the generalization of Eq. ([24])) may be written 

K n) (P cm )> = cos9\M$(P cm )) + S m0|42(fL)> (54) 

where \Mq), |Mqj), and \Xq^) are the generalizations of the states \Pq) = \Mq), 
\Pqs) = \M$l), and \P Qm3 ) = |X^°° ]) of Eqs. ©, and ©, respectively, and where 

under our assumptions that the state of the flux tube is independent of n (the adiabatic 
approximation) and that does not affect the coordinate r, 9 is independent of n. (The 
rationale for the notation (nOO) will become apparent below.) From this expression one can 
immediately obtain the generalization of our result for the elastic transition that 

e M {wf n = ^(MSp{+^)\QQ\M$\-^)) (55) 
= £g? (wr' n [cos 2 9 + szn 2 9$ M (w)] (56) 
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where ^f{w) n ' n is the valence quark model generalization of the Isgur-Wise function for 
n — > n' transitions and £,f t M (w) is exactly the same qq overlap form factor that appears in 
the elastic n = — > n' = transition. Thus the unquenched result for small w — 1 is 

p M (wf n = e Q f{w) n ' n [l - sinHpl{w - 1)] (57) 

as for the Isgur-Wise function with, once again, sin 2 6 = A qq 2 . Since the production 
of each inelastic resonant channel occurs with strength proportional to w — 1 to a positive 
integral power, the qq modification of $^ M {w) n ' n for n' > has no effect on the saturation 
of Bjorken's sum rule for £(w) to order w — 1 since it produces effects which are at least of 
order (w — l) 2 . This is in accord with the expectations outlined above that the additional 
depletion of elastic rate by qq pairs must be compensated by the explicit production of a 
(Qq) a {qd)b continuum. 

To see this we must introduce the states in the continuum orthogonal to \Mq^). To this 
end we define a complete set of states \X^-^{P cm )) in the Qqqd sector. Here n, a, and (3 
are excitation quantum numbers associated with the r, w, and v coordinates, respectively. 
These states are not the eigenstates of this sector in the absence of H^: the eigenstates are 

— * 

the |P cm ; ab(q )) defined above. However, we can expand 

\P cm , ab(q)} = £ & P \q T\X^{P cm )) . (58) 

nap 

It follows that to lowest order in 9 (or equivalently the pair creation operator) we can form 
an orthogonal set of H^-perturbed states 

\M%\P cm )) ~ \M^(P cm )) + 9Y: J d 3 q<pT\q)\Pcm;ab(q)) (59) 

ab 

\X ab (P cm ,q)) ~ \P cm -ab(q)) -eY,<pT\rT\M^(P cm )) (60) 

n 

where (a, (5) = (0,0) define the universal state of the qq pair created in a flux tube by the 
action of H^. Let us now compute the transition amplitude to the continuum: 

(Xab(+^,q)\QQ\M$\-^))^0[(+^;ri 
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- E^° 0) (? )Kl(+^)\QQ\M$l(-^))] (61) 



T / <t>T\Q)(M^(+^)\QQ\M^(-^))} . (62) 



As rriQ — > oo, P cm = mgVcm is much larger than any internal momentum so the matrix 
elements (<5|<Q<5|<5) appearing in the Qd — > Qrf and Qggrf — > Qggd transitions here are 
identical. Moreover, since the QQ current does not affect the internal state of the qq pair, 
P is required to be zero. Thus Eq. (|62|) becomes 



^ ^ n a 



(63) 

where £^ M (w) a ° is the generalization of ^9 t M (w) encountered above, namely 

^/>r° = J d 3 w ^\wye- 2m ^^f 1 } (w) (64) 

where ipi™ is the a th basis state for the expansion of the qq center of mass coordinate w. 
Expanding the exponential in powers of V cm as is appropriate for small w — 1, we obtain Eq. 
([44]) of Section II for a = and 

$ M M°* ' ^ -2zm,vL • / d 3 W ^(wTw^iw) . (65) 



Since ^/ t M (w) 00 — 1 — p^ q (w — 1) with u> — 1 = V^ m /2, to leading order in V cm the a = 

, (raOO) . 



term in Eq. ( |63"D cancels with <p^ (q) to leave 

-)) . 6 

An immediate consequence of this relation is that in the adiabatic limit (/3ft — > oo 



(X ab (+^,q )\QQ\M§\-^)} * 9 £ egf {wn% M {wr^\q ) . (66) 



£^ A/ (u>) a0 ~ /3ft 1 -> for a ^ so we have explicitly demonstrated that there is no 
nonresonant production in this limit. 
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Since we have for our discussion assumed that ip^(w) has the form of a ground state 
harmonic oscillator wave function, it is natural to use a harmonic oscillator basis as the 



orthonormal expansion functions for the variable w. Doing so, it follows that only three 



basis states give nonzero contributions to Eq. to leading order in V cm since wipf^(w) is 



proportional to the three n w = 0, £ w = 1 harmonic oscillator wave functions: 

n5/2 

^=°^ =1 > i] (w) = v^g^e-^ 2 (67) 
in a Cartesian basis, giving 

and thence 

— * — * 

)\QQ\m$\-^)) ^EC( w )T(») K= °'" =M, t K=0,4=1 '' 10) (?) • 

n,i 

(69) 

Next we note that 



where k is a positive integer by Luke's Theorem [27] so that to order V c 



(x ab (+^q)\QQ\M$\-^)) o&f ( w) ^erH K=OA=1 ' il V2 [nu ' =0 ' c=M]0) (g ) (7i) 



i 

i6\f2m q V cm 



<M<?) (72) 



Pft 

where I have introduced the notation 4> a b{q ) for the three- vector ^M"™- ^™- 1 '*] ) 

This key result has a simple interpretation. Nonresonant production at order w — 1 
requires that the current QQ acting on Qqqd create neither the ground state nor the P-wave 
resonances. However, it cannot excite the qq internal coordinate v and, if it excites r then 
to order w — 1 it has just produced the Qqqd component of either the ground state or the 
P-wave resonances. Hence nonresonant production to order w — 1 occurs purely by excitation 
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of the qq coordinate w to £ w = 1. The factors 4> a b{(f ) are simply the projections of these 
£ w = 1 states onto the continuum eigenstates consisting of mesons a and b with relative 
momentum q. 

We are now in a position to verify that the Ap^ ea contribution to the slope of the Isgur- 
Wise function is indeed compensated by the production of these continuum (Qq) a (qd)b states. 
The probability for their production at w is, up to (w — l) 2 corrections, 

dP{P Q -> continuum) ~ £ j d 3 q\(X ab (+^, q)\QQ\M§\-^))\ 2 (73) 

ab " 

which since £ ab / d 3 q<j)i b (q )(jP ah {q ) = 6 ij , gives 



. 26 2 m 2 a V2 , . 

dP(Pq — > continuum) ~ -| (74) 

(75) 

Pft 

for the flavor q. On the other hand, according to Eq. (|53D, the contribution of flavor q to 
the loss of elastic rate is 

Am? 

dP(P Q - P Q + V Q ) Q - qqd ~ --^^ 2 (^ - 1) • (76) 

The two rates match, explicitly demonstrating the connection of Ap 2 ea to the nonresonant 
continuum. 

3. Production of Exclusive Nonresonant States at Low w — 1 

It remains to assess the fractional population of individual continuum channels inside of 
the total given by Eq. ((flf). To do this we must calculate (j) a b(Q )> which from its definition 
is 

4>%*>(?)P(P' m -P m ) = (P' m ;ab^)\X^(P m )) (77) 

for the case {na(3) = (0[n w = 0,£ w = 1,*]0). These calculations are straightforward, but 
would be quite tedious without the introduction of several tricks described in Appendix 
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C. Results of the calculations of T{% (w) and t^ 2 (w) for number of low-lying nonresonant 
channels are given in Table I. 



III. UNQUENCHING HEAVY QUARK DECAY: RESULTS 

We now turn to the quantitative evaluation of A-p 2 sea (which is a reflection of total non- 
resonant production) and then to the distribution of these decays into exclusive nonresonant 
channels. 

A. The Total Nonresonant Rate 

As shown in Appendix B, the light quark amplitudes r] ua = rj^ can be determined from 
strong decays to be 

rj qq - ~ 0.9 GeV (78) 

with an uncertainty of a factor of two mainly arising from a strong dependence on the poorly 
known quantity (3 pc . It follows from Eq. fl35|) that 



c qq ~ 0.5 . (79) 



Assuming SU(3) symmetry for the contribution to Ap^ ea , one then obtains from Eq. fl53|) 

ApL * \ , (80) 
corresponding to an increase of p 2 from the value 0.74 ± 0.05 quoted in ISGW2 to a value 



near unity. Either value would be in reasonable agreement with measurements jl5]. Via 
the Bjorken sum rule, such a p 2 sea would be consistent with the possibility discussed in the 
Introduction that 16 ± 8% of the inclusive semileptonic B rate is in nonresonant channels. 

This reasonable quantitative correspondence between our calculated Ap^ ea and a pos- 
sible experimental anomaly should not be taken too seriously. The missing non-(D + D*) 
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rate attributed to nonresonant production might be due to an ISGW2 underestimate of 
excited resonance production ||18|| , or to an experimental overestimate of non-(D + D*) pro- 
duction. Moreover, while it is a "canonical estimate", Ap^ ea is subject to very substantial 
uncertainties: see Table II. 

B. The Rate to Low-Lying Exclusive Nonresonant Channels 

Using the amplitudes of Table I and the formulas of Appendix C, one can easily calcu- 
late the fractions of Ap^ ea due to the individual low-mass nonresonant channels shown in 
Table III. Note that the nonresonant rate is highly fragmented: none of the many channels 
tabulated account for more than 4% of the inclusive nonresonant rate (and thus no more 
than about 1% of the total semileptonic rate). These results are consistent with previous 
studies of single low energy pion emission using heavy quark chiral perturbation theory ||11|| . 
Note also that the thirty final states considered here account for only about 40% of the 
nonresonant rate, so that most of this rate resides in states ab composed of highly excited 
resonances. 

There are several reasons why the results given in Table III must be interpreted with 
caution. The most prominent is simply that they depend roughly on /3~ c 5 ! Recalling that 
P pc = n l / b (3 pc , our canonical value for (3 pc is 0.4 GeV, but over the range of (3 pc allowed by 
Table II, it varies from its canonical value by ±0.1 GeV. Although the factor n 1//5 makes 
this range narrower than that of (3 pc itself, it still leaves us with an uncertainty of more than 
a factor of two on the basic unit for production of exclusive nonresonant channels. 

The second word of caution concerns large corrections to the itlq — > oo limit studied here 
which arise for real B decay. The problem is phase space: in the heavy quark limit, all final 
states occur at the w of the underlying quark production process, i.e., the full mass range of 
the final state spectrum of states is negligible compared to the heavy quark energies even at 
low w — 1. In actual B — > X c ivi decays, each final state X c will have a Dalitz plot that is a 
shrinking fraction of the b — > civt Dalitz plot as mj c — > mg. Thus in the heavy quark limit 
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the loss of rate from the elastic channels D + D* would be locally compensated in the variable 
w. For finite m b and m c , however, a loss of elastic rate will still occur via Ap 2 ea , but the 
compensating channels will experience a delayed turn-on because of their thresholds; indeed, 
some processes which would have helped to compensate Ap 2 sea will be kinematically forbid- 
den. These phase space suppressions lead to Agc^/mQ-type corrections to the inclusive 
rate, and therefore also corrections to the accuracy of quark-hadron duality. However, it is 
known on general grounds from heavy quark symmetry and the operator product expansion 
(OPE) that, as the energy release m;, — m c — > oo, the leading corrections to the inclusive 
rate must be of order K^ cd /itlq |28|j29|1 . The resolution of this apparent paradox has been 



discussed in Ref. [[Kj: the OPE has a radius of convergence which does not include the region 
in which (significant) hadronic thresholds are turning on, so for real b — > c transitions there 
can be Agc^/mg corrections, and they can be substantial. However, associated with these 
Aqcd/^q threshold effects, which would diminish the integrated contribution of individual 
hadronic channels below the level required for perfect duality, are A.QCD/ m Q corrections to 
the rate for the production of such channels which enhance their production once they are 
above threshold. These counterbalancing effects soften the breaking of duality: they are the 
precursors of the perfect cancellation of Ag C D/ m Q effects that must occur as m b — m c — > oo. 
How then should one make use of Table III for real b — > c transitions? The amplitudes 
and 7-3/2 shown are the leading order predictions for their respective channels. For the 
reasons just outlined, we can expect that they will be enhanced by Agc^/mQ corrections 
which are expected to be of the order of 25 — 50% [|30|. However, despite this effect, I believe 



that the dramatically reduced size of their Dalitz plots (relative to that for the underlying 
quark process b — > due) will reduce the actual population of all nonresonant states well 
below that expected from Ap 2 ea . Indeed, each channel shown in Table III has a continuous 
spectrum of masses from its threshold up to masses exceeding tub- Consider, for example, 



the D+p channel. Its mass is given by m£, P (q 2 ) = \Jm 2 D + q 2 + ^Jm 2 + q 2 , while its dominant 
.D-wave production rate is proportional to q 4 e~ q ^ if3 , where 2/3 ~ 1 GeV. The production is 
therefore very weak at low masses where the available phase space is generous, and peaks 
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at mop ~ tub where it vanishes. Given these basic kinematic facts, it is clear that even 
this simple exclusive channel will be produced at a rate far less than in the heavy quark 
limit, and that most of Ap^ ea will be uncompensated. The third column of Table III gives 
the phase space factors by which individual channels will be reduced in real B — > XJLvt 
decays relative to the heavy quark limit; the fourth column gives a very rough estimate for 
the net suppression of each channel as a product of this phase space factor and a generous 
guess that, after KqcdIthq corrections, each channel has a compensatory increase of 50%. 
Considering that in general the untabulated channels of yet more highly excited states ab 
will suffer even greater phase space suppressions, an overall diminution of the nonresonant 
rate by at least a factor of four seems likely. 



IV. CONCLUSIONS 



Although the successes of the valence quark model and the arguments of the large N c 
limit provide indications that sea quarks play a relatively minor role in hadronic physics, 
this hope is far from being justified by our current understanding. Some failures of the quark 
model (e.g., the proton spin crisis) and the known existence of strong real and virtual decay 
channel couplings indeed make blithely ignoring the role of qq pairs both phenomenologically 
and theoretically untenable. In this work I have examined the influence of qq pairs on the 
simplest "real" hadrons: heavy quark mesons like the B. 

This study has led to a number of qualitative insights which I believe are quite general in 
nature. In earlier work on "unquenching the quark model" , the success of the valence quark 
model in spectroscopy was shown to have a possible basis in the validity of an adiabatic 
approximation. In this approximation, both the confining flux tube and the many qq pairs 
it generates remain in their adiabatically evolving ground state as the valence quarks move. 
In this work I have shown that the same approximation leads to valence quark and therefore 
resonance dominance of the simplest current matrix elements: matrix elements of 

heavy quark mesons. The physical picture behind our results is simple and appealing. 

27 



According to heavy quark symmetry, at small recoil the Qi — > Q2 decay of the Qi ground 
state Pq 1 will lead with unit probability to Pq 2 and Vq 2 no matter how complicated the 
QCD "brown muck" might be. This simple observation makes it clear that for Q2^Qi matrix 
elements the issue is not the probability of qq pairs in Pq 1 but rather how rapidly asw-1 
increases these pairs fail to overlap with those in Pq 2 , Vq 2 , and the Q20I excited states. I have 
shown explicitly that in the adiabatic limit these overlaps are perfect so that only valence 
states (the resonances) are produced. Moreover, I showed that violations of the adiabatic 
approximation can be directly associated with the production of nonresonant states. Thus 
this study leads to a way of understanding how the valence quark model can be so successful 
even though hadrons are full of qq pairs. 

While they are quantitatively very crude, these calculations also have interesting con- 
sequences for real B semileptonic decays. First of all, they suggest that p 2 dyn = p 2 — \ is 
composed of two comparable parts: p 2 resonant ~ \ and p 2 nonresonant ^ \, though we have 
stressed that this split of p 2 dyn is very model-dependent. In the heavy quark limit, Bjorken's 
sum rule would then lead one to expect (using the central experimental value of the (D + D*) 
fraction) that roughly 24% (with a 50% error) of semileptonic decays go into resonances (both 
ordinary Qd mesons and Qd hybrids), and roughly 12% (with an error of a factor of two) 
go into nonresonant states. Since the former states will have most of their strength in the 
2.4 — 3.0 GeV region, they will suffer, relative to the heavy quark limit, phase space suppres- 
sion factors varying from only 0.75 to 0.50 over this range which may be fully compensated 
by the ^QCD/ m Q enhancements required by asymptotic duality p9|p0| . In contrast, we have 



seen that nonresonant states are expected to populate very high masses peaking in strength 
near m# and so to suffer a substantial net suppression factor of at least |. From this study 
I therefore expect that ~5% of B semileptonic decays will be nonresonant! 

As a corollary of this last observation, I note that if a 12% nonresonant semileptonic frac- 
tion is required for duality but only a quarter of this is realized, then duality will fail from 
this effect alone by ~ 10%, as anticipated in Ref. |30[|. There is, however, a minor inconsis- 



tency associated with this conclusion. Experiment requires that 36 ± 6% of B semileptonic 
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decays be non-(D + D*) decays, in contrast to the ~ 25% we would have estimated from 
the preceeding. As mentioned earlier, this could simply mean that the ISGW2 model under- 
predicts the production of excited charm mesons QXHJ or that experiment has overestimated 
non-(.D + D*) production. Determining whether this discrepancy is real will require a more 
quantitative calculation than this one (and probably additional experimental measurements 
as well). 

Detailed experimental studies of the structure of the hadronic final state in semileptonic 
b — > c decays can therefore answer some fundamental questions about the role of qq pairs and 
about duality in strong QCD. A vital feature of these systems is that duality is underwritten 
by Bjorken's sum rule, requiring an exact and local duality between quark- and hadronic- 
level decays in the heavy quark limit. In particular, the experimental determination of 
the strength and structure of these nonresonant contributions would immediately test the 
conclusions reached here that these qq effects are highly suppressed in real b — > c decays, that 
such decays extend to very high masses, and that they are highly fragmented into many small 
channels. Independent of the outcome, examining this problem in Nature's simplest hadronic 
system under the action of its simplest current (a heavy-to-heavy nonsinglet transition) 
should prove to be an excellent starting point for eventually understanding the qq sea in 
all strongly interacting matter. In particular, given the complexity of QCD, this seems an 
essential first step before tackling the problems of duality and nonresonant production in 
ordinary deep inelastic scattering. 
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APPENDIX A: Flux Tubes and A Critique of Nonrelativistic Kinematics 

It is not difficult to make the simplification x/jf t (w,f) = ipf t (w) of Eq. ( pT|) more ffirx- 
tube-like. In the case where ipf t depends on w and r, Eqs. ((43])- ([44]) become 

— * — * 

d 3 wd 3 r^j* QlI {f )V»/ t K f ) e - 2im ^-^e- imdV?cm - r > Q( j(r )ip ft (w, f ) (82) 



or, defining 



we have 



md = ^J d 3 wd 3 r(2m q w + m d r) 2 \tp Q j(r)-il) ft (w,r)\ 2 (84) 
= ^((2m q w + m d f) 2 ) (85) 

which reduces to the simplified results of the text in the appropriate limits. 

Now consider a generic example of a more realistic ipf t that has a flux tube's shape: 

M^r) = &L e -W*<t(iB-r) (86) 

7T ' 

where f = r/r, w± = w — (w-f)f, and t(w-r), which depends only on the longitudinal variable 
w ■ r, is a normalized tube-like function. (E.g., one might have t = -^0(w ■ f)6{r — w - f) to 
create a cylindrical wavefunction that is gaussian transverse to f and constant between Q 
and d) With such a wavefunction 



The first term is as expected intuitively: it is the unchanged w± part of the result of the text. 
One might also naively interpret the second term as the old r term plus a new longitudinal 
contribution due to the assumed spatial distribution of w ■ r in a tube-like configuration 
along r. 
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I believe that the physics is more subtle than this. Consider the origin of p^y in the 
nonrelativistic kinematics of our model. In the Qd sector (see Eq. ([]])), fq = —mdf/mq 
in the center of mass frame as mq — > oo, while in Qqqd (see Eq. (0)), fq — ► — (2m q w + 
mdX )/rriQ. Since nonrelativistically p 2 = |mg(rg), we see that Eq. (|50| ) for p^j and Eq. 
(|S5|) for pq qq g are simply consequences of these nonrelativistic relations. 

To see the dangers of this approximation when the string tension and its renormalization 
are large compared to m q , consider the calculation of the mass of a system of heavy quarks 
Q and d at separation f connected by a renormalized flux tube, i.e., of the state (p4|) of the 
text which has string tension b^ f since it has the appropriate admixture of qq pairs. If its 
mass were determined nonrelativistically one would obtain 

M nr -m Q = m d + — *M (88) 

1 + C <7<? 

i.e., the effective mass opposite Q (against which it must recoil to conserve the position of 
the center of mass) would be the probability-weighted masses of the pure d state and the 
qqd admixture. On adding interactions (both the diagonal potential b$r and the off-diagonal 
perturbation which mixes Qd and Qqqd), we obtain the correct answer 

M adiabatic _ mQ = ^ + ^ _ (gg) 

Thus the mass 2m q does not in this circumstance have an independent reality as indicated by 
the nonrelativistic kinematics just described: it is subsumed into the properly renormalized 
string tension. Indeed, it is an implicit assumption of the model for the qq pairs described 
in the text that the unquenched flux tube also behaves like a relativistic string: it should 
support (only) transverse waves moving with the speed of light. Thus while the qq pairs 
change the longitudinal distribution of fq, this effect is already described in the flux tube 
model by the mass b^ f r residing in the flux tube: when r increases, not only does d move, 
but so does the center of mass of the flux tube. It would therefore be double-counting 
to include the effect on (r%) of the 2m q w ■ f term of Eq. fl57p. I should hasten to add 
that the nonrelativistic quark model used in this paper does not explicitly take this effect 
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into account. To do so would require a full treatment of the flux tube degrees of freedom 
(versus the adiabatic, potential-model approximation used here). Nevertheless, even the 
nonrelativistic quark model has undoubtedly already taken some of this effect into account 
implicitly by its choice of such free parameters as ma- {E.g., in many applications the 
constituent quark mass is effectively ro^ = mJJ + ~&(r)). In contrast, there is no mechanism 
in the quark potential model to take into account transverse displacements of Q relative to f 
|3j| . These transverse displacements are the true degrees of freedom of the (quenched and 



unquenched) flux tubes and the reaction of Q to them makes a non-potential-model-type 
transverse contribution to p 2 . 

By renormalizing the string tension from 6 - ► ^>N f < &o> QQ pairs increase the longitudinal 
contribution to p 2 (at least in the nonrelativistic approximation << m^). However, this 
increase is not compensated by nonresonant production, since it is this same string tension 
6jv> < b which controls the structure and thereby the production of excited resonances. 
I.e., the longitudinal effect of the pairs is real, but it simply renormalizes resonance physics. 
The dynamics behind this balancing act, characteristic of the adiabatic approximation, can 
be seen by calculating the contribution of the flux tube to the Isgur-Wise function |TT 



since the flux tube has only transverse internal degrees of freedom, it has no impact on 
longitudinal overlap integrals over the total separation r. Based on the preceding arguments, 
the transverse contributions to (tq) would also be those of a relativistic string with string 
tension b or b^ f in the quenched and unquenched flux tubes, respectively. In contrast to the 
longitudinal contribution to (tq), the transverse contributions correspond to the reaction of 
Q to real internal degrees of freedom, and these degrees of freedom (both gluonic and qq) 
can be excited by the action of the QTQ current. When acting on the pure Qd piece of the 
state, the current excites hybrid mesons which in the quenched limit exactly compensate for 
the loss of rate from the elastic channel due to transverse contributions to (tq) In the 



Qqqd sector, the current could in principle excite either of the strings internal to mesons 
{Qo)a or (qd)b, or it could excite the center of mass of the qq pair. In the quark potential 
model approximation to this latter process, which is of course the one of interest for this 
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paper, one would recover the result shown in Eq. ( j37D less the longitudinal part of w: 



responding to two of the three degrees of freedom of w being active. Note that for this 
picture to be consistent, the total transverse contribution to (tq) must be that of a relativis- 
tic string with string tension b^,; the decomposition into Qd and Qqqd components is only 
useful in identifying the compensating channels required by Bjorken's sum rule. However, 
in the renormalized string picture (rq) of course depends just on , while Eq. ([)]]) shows 
a contribution proportional to 2m q . It would be interesting to examine how the dynamics 
of pair creation in the flux tube leads to such a term. I speculate that the mechanism is 
the "consumption" of a piece of flux tube of length Ar ~ 2m q /b in a nonlocal pair creation 
process. 

In summary, in this Appendix I have described several subtleties in the description of 
qq pair creation in the flux tube model, and pointed out some interesting issues which arise 
in the physics of the renormalized flux tube. While I believe these matters are important 
conceptually and are worthy of further study, I am convinced that other uncertainties de- 
scribed in the text are of far greater impact numerically on our results. Given this and the 
great convenience of the spherical approximation, I therefore chose this simpler if less basic 
framework for the discussion of the text. 




The first term of this formula differs from the expression of the text by the factor | cor- 
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APPENDIX B: Determining r] qq 

Eq. ( pi) defines the action of the pair creation Hamiltonian on I-Pqj). This perturbation 



not only produces pairs to make the eigenstate \Pq) of Eq. fl24|), but also leads to strong 
decays. In particular, the projection of the state (^) onto the continuum states \P cm ; ab(q )) 
determines the Pq — > (Qq) a (<ld)b coupling constants. By heavy quark symmetry ||, the 
same dynamics determine the Pq — > (Qq) a (qd)b coupling constants, where Pq is the vector 
partner of the pseudoscalar state Pq. In this Appendix I use these facts to determine the 
strength parameter r] qq - of Eq. ( |17D by comparing to the decays Pq — > Pqti. 

I begin with a practical matter. As itlq — > oo, the decays Pq — > Pqtt are forbidden 
since Pq and Pq become degenerate heavy quark spin partners. However, T(K* — > Kir) is 
known, and T(D* — > Dtt) can be deduced if one assumes that the successful phenomenology 
of magnetic dipole decays can be extended to T(D* — > -D7). (The total width of the D* is so 
small that only decay branching fractions and not decay widths are known. If one takes the 
K* —>■ Ktc and K* —>■ A7 decays and scales them appropriately in mg assuming that heavy 
quark scaling works all the way down to m s , the observed D* — > D7r and D* — > branching 
ratios are explained nearly perfectly. This is another example of the often-noted fact that in 
many circumstances a strange quark behaves like heavy quark.) Since the branching ratio 
for D*° -> D°n° is well-determined experimentally, I will use the value T(D*° — > D°ir ) ~ 30 
keV deduced in this manner as input "data" . 

The calculations themselves are simple. If g^D^d^D *^ is the effective Lagrangian 
density for the decay leading to T(D*° — > D°tt°) = 48 9 n ^2 , then in the center of mass, with 
the pion emitted with momentum q from a D* with polarization +1 along z, 

im± fj m^fofPfiffi r r r 

7= 7- = Vos\ a — / d r / a v / a w 

v / 2(2vr) 9 / 2 qq \ 3 8vr 6 J J J 

c -y^+i^-yi(r-a+i^-i^-h^ ^ 1 ^ c -iUf+a- y)^ pc (if) + (92) 

in which the for the pure n° decay via r] m has been explicitly included, fax is the 
"mock meson" mass given by the sum of the constituent quark masses, and ip pc (v) + = 
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v + e iPpcV 2 . The integrals are straightforward and give 



9o = -ifar — To\ izo* — " UU> e^ (93) 



where 



Wo'Pl ' + Pi i lhTW ) + p! { 2- a - b) (94) 

Add = + a 2 (/% + $ + $ t ) + (6 2 + !)(/% + ^) + a(3 2 D + (b - 2ab - a)[3 2 (95) 

C = + # (96) 



= + ^ + S%r 2 (97) 



2[(/% + + # + (98) 

bDD 2[(fi, + + % + ft] ■ ^ 

From the "measured" D°* — > D°7r° width, we can deduce that g ~ 11. With our 
canonical parameters it follows that, with f3 pc expressed in GeV, 

^-^0.32[^] 3 / 2 [4l±^ ]5 / 2 (100) 

where since f3 D ~ (3 n ~ (3ft/ f 1 ^ 2 I have been able to explicitly display the dependence on / 
as well as j3 pc . We see that for a variation of ±1 around the canonical value / = 2, r) qq - varies 
by less than 10%. Thus we conclude that (since f3 2 c » 0.06 GeV 2 ), 

"-°' 9 loifevl V2GeV (101) 

i.e., r/gg ~ 0.9 for the canonical value (3 pc = 0.58 GeV. We also see from this formula the 
expected result that as the pair creation operator becomes more pointlike, r) qq - — > oo. 
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APPENDIX C: Calculating t$ and r 3 ( ^ for Selected Low- lying Exclusive Non- 
resonant Channels 



As discussed in Ref. |[17| , the semileptonic decays B — > D^ 1 }, £ug are governed in the heavy 
quark limit by generalized Isgur-Wise functions which determine all of the form factors for 
the decay of the B with s^ e = \ 1° both of the states of a heavy quark spin multiplet with 
quantum numbers s'\ l . As described in the text and elsewhere [|17]], 

p 2 = - + Ap 2 pert + (? Qi + Ap 2 sea (102) 

where the \ is Bjorken's relativistic correction fTS] , Ap 2 ert is a perturbative QCD radiative 
correction, and p^j and Ap 2 sea are the contributions to the slope of from the valence 
and sea quarks, respectively As we have seen, these latter two contributions may be related 
to the rates of decay into inelastic channels by 

p 2 q ,= E k$| 2 + 2 E l^l 2 (103) 

1 _|_ - o _|_ - 

m= 2 Qd resonances p= ^ resonances 

Pla= + E |T^| a + 2 + E l^l 2 (104) 

m=i continuum P = f continuum 

where the r's are the appropriate Isgur-Wise functions. Once the r's are specified, all 
transition form factors to the states of a heavy quark spin multiplet may be determined from 
symmetry considerations. For this reason, it is useful to calculate the r's in the simplest 
possible setting, namely for the case where b and c are spinless. This is possible since the 

r's depend only on the dynamics of the light degrees of freedom, i.e., on the transition 

Q *t _ i- _. j-k'i — i+ nr 3 + 
b i — 2 h e ~ 2 UI 2 • 

Another key simplification arises from the dynamics of the pair creation process. As 
demonstrated in the text, when b(v) — > c(V) in a Qqqd state, excitation of the variable r 
cannot lead to a contribution of order (w — 1). Furthermore, the variable v cannot be excited 
since this is a qq internal coordinate. Thus to contribute at order (w — 1), b(v) — > c(v') must 
kick the w coordinate into an £ w = 1 state: such a state is the parent of all nonresonant 
production to order (w — 1). Thus Ap 2 sea arises entirely from the "decay" of the lowest 

36 



£ w — 1 excited state of cqqd arising from the b — > c transition from frggd. With the pair 
in J p = + , the decay can thus occur from the six states w + j, w + |, u> T> w o I) w_ T, 
and u>_ I, depending on which component of w is excited by the recoil. Here w + ,wo,W- 
represent the components of the £ w — 1 state, and j, j represent the spin state of the d 
spectator quark. Since the total decay rate of bqqd j and frggd j must be the same, we can 
simplify if we average rates over the initial d spin and over directions of P cm (or equivalently, 
over directions of w). Then since the average over the six states just listed will be the same 
as the average over the two j = ^ and four j = | states formed from them, we can 
deal with "parent" states that are states with good angular momenta and which therefore 
uniquely feed the s'J 1 = \ + and | + states, respectively. 

Let me provide an example: the production of the (D + D*)n nonresonant states. Since 
we know that all of Apj. ea arises from the "parent" state, the fraction of Apj. ea coming from a 
given channel can just be obtained as the jrrij average of the square of an overlap between a 
given jrrij in Qqqd and the two particle continuum state of interest. Thus from the jrrij = \\ 
state we can extract ((D + D*)i-iir(q )\cqqd] | + |)' which ought to leave the (D + D*)n 
system in an S-wave. Explicitly, as ttiq — > oo, 

A. 1 + _ f j3... f * f jzPT^-hPAw+ht) 2 ^ 



((D + D*) h -Mq)\cqqd-- -> = J d*w J d 3 v J dV-^e^ 



7,3/4 



1 „-iUr+w-\v)P% 2 „-k0Lv 2 ^ft 2 -1/3 2 h w 2 Pb 2 -i/3%r 2 



(27r) 3 / 2 7T 3 / 4 7T 3 / 4 7T 3 / 4 

H(w,v) (105) 



where 



s = (T \J\(U - iT)iy|(Tt «- - [Tl + ITK- II v + )^(-w + I -«;, T)) (106) 

is the spin overlap matrix element of the three quarks qqd, respectively. We get 

((D + D*)i-in(q )\cqqd; \ + \) = -7°°- - 4°° (107) 

2 2 Z 2 

where 
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( r lv I ( P r I d3 r J^ e -*§-(^-^-^ J 



(108) 



with 



E 2 = (3l(w + \vf + ft{f- w + \vf + /?> 2 + /3> 2 + (3^ ■ 



(109) 



In these formulas I have distinguished between /3f> and to allow the "ancestry" of terms 
to be traced, even though (3b = (3d from heavy quark symmetry. This integral is easily 
done, giving 



rOO 



rOO 



c 00 5- + c u 



ft 



-g 2 /8/3 



(110) 



where I , d 



oo „oo ^oo 



and (3 are given below. For the problem at hand, we get 



((D + ^*)i-i7r(g )\cqqd; \- h = - 

2 2 Z Z 



3c 00 + c! 



.00?+ , ^00? 



lJ A 2 

2 



y ^2 



o 00 , 00 S_ 

j PI 



o 00 , 00 Q_ 

3c +C(, 



A 

j00 e -g 2 /8/3 
2 



j00 e -g 2 /8^ 



J00 e -V/ 



^00 (fig) 



111) 



a pure S'-wave decay as required, with partial wave amplitude Ay 2 



3c 00 + 



joo e g 2 /s/3 _ Note that with Y 00 factored out, |v4i/ 2 1 2 is already the probability for 



this channel integrated over angles Q q , leaving only an integral / dqq 2 \Ai/2\ 2 to be done to 



sum over all (D + D*)ix states with quantum numbers at any fixed value of (w — 1). 

Next consider ((D + D*) 1 - i 7r(g )\cqqd; f + f), which proceeds by replacing \/f(— w + | 

22 ' " 



I)) by — w + t in Eq. (|106|) . With this change one gets 



22 z z 



'_ T 00 

2 z+ 



3 Q() gzg+ TOO -g3/8^ 



^[V4^ c °°|,j 00 e -^/^]y 21 (fi 



9/ ' 



(112) 



a pure D-w&ve decay as required. Moreover, since —J\ is the Clebsch-Gordan coefficient 



for coupling (D + D*)i-i and an £ = 2, m = 1 pion into a M state, we can deduce that 

22 
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decays to this whole angular momentum multiplet with (D + D*)tt in a D-wave coupled to 
s 't 1 = I are controlled by a D-wave amplitude A 3 / 2 = — \/A^cf-^I 0Q e^ q / 8/3 . 

To complete this pedagogical example, I note that since the partial wave decay amplitudes 
are independent of the magnetic substate m, 

\i:\aim 2 = 1\a 1/2 \ 2 (us) 

\T.M = \M , (114) 

where the | arises from averaging over the six jrrij states. On comparison with Eq. (|104|) , 
we see that t 1/2 = ^A 1/2 ^lAp 2 sea and r 3/2 = ^A 3/2 ^Ap 2 sea . Note that \Ap 2 sea appears 
since the overlap amplitudes Ay 2 an d ^3/2 as calculated are for a single flavor, while the 
factor y| arises as a residue of the angular and spin averaging. 

The tricks outlined here are more powerful for more complex decays. I will give one 
partial illustration: ([(D+D*)p\ 3 + 3 \cqqd; | + |), where the subscripts on the bracket are total 
spin quantum numbers, but do not include relative orbital angular momentum. The overlap 
integral for this decay is obtained by replacing the pion spin wavefunction ^/|(ti — IT) by 
TT, \J\{-w+ I -w z t) by -w + |, and ^ by (3 P in Eq. fll06D to give 

22 2 2 

= _ v ^[ 2c 00 + c 00^] / 00 e - g 2 /8r 

Pv 

= _y|vS([3c« + c«j]r m 

By examining the Clebsch-Gordan coefficients for coupling the spin state [(D + D*)p]z + 3 

2 2 

to a relative orbital S"-wave or a D-wave to get a || state, we can deduce that the 
S- and D-w&ve amplitudes for this decay are — y^[\/47r(3c 00 + c°°|^)/ 00 e"' ?2 / 8 ^ ] and 

2 _2 

[v47rc^ %nl m e~ q 1^ ], respectively. One very simple overlap integral thus gives two par- 
tial wave amplitudes simultaneously. 
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A complete set of results are given in the text in Table II in terms of the following basic 
integrals: 

^oo Q&j ' 



00 



J 00 = / 

«7 



rlO _ rlO 
A ijk - A 



PI 



V/8/3 



10 QiQjQk 



c t qi5 jk + Cj qjSik + c k q k 5ij + c ijk ^ 



jOl 

ijk 



r01 



cfqibjk + cfqjdik + c^q^j + c£ fc ^ 



,01 



01. 



2 

oi qiqjqk 



yv a 



where J™ is specified in terms of 



r 00 



5/2, 



= +4aA 



c 00 



a 

v 2 



/3l + 2^ o/% 1 



&) 5 



& 2 = + « 2 (/% + Pi + P%) + & 2 (/?l + Pi) 

+\{Pl + Pi) + aP 2 D + (b- 2ab - a)Pl 

2 i o2 



P 



J w 



Pd + P 



ft 



PlPl 
Pl + Pl 



(P 2 B+Pl)(P 2 D-Pl)+^ 

" 2 [(Pb + PV){Pd + PZ + P%) - Pi] 

, P 2 M + P%) 

' 2[(/3| + Pl){Pl + PI + P%) - Pi] 

l = l + J.( /% + 2/£ )2 j_ ( i_ a _ 
Similarly, I}® k is specified in terms of 

^pTpi I2 p% 2 ~p% 2 



bf 



I w = 



40 



40 



r io _ 



Z 5 



StT 3 / 4 /?^^^ P w ^ Py^ 



,P 



2 
it;** 



2 
if** 



0^ + 201 1 1_ 

Pl + Pl 2 7V 2 

(1 / 1 \ / 1 , /^', 

P'l 

1 w /%+jg° 



1 1 



^** ) 



2 



~ 2 
Pv** 



(116) 
(117) 
(118) 



(119) 
(120) 
(121) 

(122) 
(123) 
(124) 

(125) 
(126) 
(127) 



(128) 
(129) 
(130) 
(131) 
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r 10 = +- 



,1 



o** ) 



" 2 3 2 

ti** /A^_B 



2/33. 



Pw** Pb + ^ 



-(- + 2a**)(- - a** ^**) ( ^,9 zoo ) 
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^2 



1 1 

+a**(- + 



2 
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3 • ^/:* 



~ 2 
Pv** 



(132) 



where /?„**, /?«,**, a**, &**, and ^- are given by the formulas for the variables /3 2 , 



fty, /3 2 , a, b, and -p, respectively, with p D 
Finally, is specified in terms of 



(3d** everywhere. 



rOl 



.01 



.01 



.01 



.01 



^pTpTp^pt^ 2 
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^ 2 



a a /3 2 ,/£ + 2/£ w l 
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y p 2 B +py K 2 
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where /3„ a , $, 2 , /?^ 2 , Oa; and J^- are given by the formulas for the 1^9 variables /3 2 , /3 2 , 

Pa 

a, 6, and -p, respectively, with p n — > /3 a everywhere. 



T^v^ -a„ + 6 a )(- 

2^2 



(133) 
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TABLE I. Ty^ and t^) 2 f° r nonresonant semileptonic decays to low-lying exclusive chan- 



nels in units of as = 

(3™ = V4^(4 g)/ 10 e-^ 2 / 8 ^ 



4¥(3c 00 + c^q 2 /^)I 00 e-i 2 / 8 ^, 



OLD 



I^(c^q 2 /^)I 00 e-i 2 / s ^, 



10 



and 7° 1 = v / 47r(c? 1 fc g 3 / /3^ a )I 01 e~ q2 '^ s ^ a , as denned in Appendix C. Shown explicitly are the ampli- 
tudes for emission of a (D + D*), D^* 2 , or D** 2 and a positively charged light hadron (ir + , p + , a 2 , 
af, b±, or Oq) from B° by uu pair creation. (Here (D + D*), D^* 2 , and -D** 2 are the lowest-lying 



jK _i-3- 



2 , | , and g heavy quark spin multiplets, respectively.) The subscripts on a channel in 
the first column are the total spin (the s^ a of the charmed meson a plus the spin of the light me- 
son b), and the ab relative orbital angular momentum, respectively; the second and third columns 
define the total s 1 ^ 1 as either ^ + or | + . Note that since partial wave amplitudes with respect to 
the direction of the vector q are given, full rates to a channel are obtained by integrating over q 2 dq 
and not d 3 q, i.e., a factor y/AHi is included in each amplitude. Shown in parentheses under each 
allowed amplitude is the fraction of the uu rate going into this channel in percent, based on Eqs. 



Hap and (114) 



channel 



(fraction of uu rate (%)) 



(channel) 
3/2 I 

(fraction of uu rate (%)) 



[(£> + £)>+]!. 



[{D + D*)ir+]x D 



l(D + D*)p+] 



-a s 
(0.4) 



-an 



(1-5) 



4 

3 a S 



(1.0) 



[(D + 



[(D + D*)p+] hs 
[(D + D*)p+h D 



(2.0) 
(0.1) 



(2.0) 



+ \/k a D 



46 



[- D 3/2 7r+ ]fP 

[£*/>+] fp 

l D t/2P + kF 
[DT/2P + h 2 F 

[DT/2P%P 
[(D + D*)a+] i p 

[(Z> + Z>>+]sp 



(0.5) 



+v /§ - 3/3/° - 2/3H - ^ 



(1.1) 



/ 16 10 

5 ' 



(0.1) 
(1.0) 



(0.3) 

+^?[/?)°-A 10 ] 
(0.3) 



+ V / |[5/5] + /3/ + 9/3n + ^27 10 
(0.8) 

+Vf[3/3f - A 01 + Wf ] + Vf 7 01 
(2.4) 

+ +A° 1 ] + V / |7 01 

(1.1) 



(0.5) 

V 3^' V 15 ' 

(0.8) 

+v / ^ 10 

(0.1) 

-vf^f + s/^H-vfT 10 

(0.4) 

+yT2[/3j°+/3n + \/i7 10 
(2.6) 

32 10 
25 ' 

(0.1) 
(0.1) 

+vlH 10 

(0.1) 

-vl^) - #°1 

(0.1) 



9 y^3 
(1.3) 


(0.0) 

+v /|[2/3jo - 4/3/0 - 9$°] + vV 
(5.6) 

+ V / |[2^ + 3/3n + V / |7 01 
(1.0) 

+ v /|[ / 3f - 2/3f ] - y^/ 1 
(0.3) 
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[{D + D*)ath_ F - +^h m 

(0.0) 

\{D + D*)a+] i p + V / |[/3f + /3f + 3/Jg 1 ] + v^ 01 + V / |[2/?f + + 3/3f ] + yfi 7 01 

(2.1) (2.1) 



3 l^a^j i "^i J i y 3 I V 3 V 15 

(0.6) (0.7) 
[(D + D*)bt]s F - +47 01 

(0.1) 

[(^ + D*)bt]i P +y /l\ff - 3/3^ - Wf] ~ Vf 7 01 "Vf [2/?f + Wf] ~ Vf 7 01 

(0.3) (3.3) 
[(D + D*)ath 2P - +yi2[/?f +A 01 ] + V / i7° 1 

(2.5) 

(0.1) (0.1) 
l(D + D*)a+]s p -y/ff? - ^ +yi[3/3f - 2/3f ] + yfei 

(0.0) (0.6) 



01 



(0.0) 
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TABLE II. Assessing the Uncertainties in c qq - and Ap 2 ei 





Pft (GeV) 


(3 pc (GeV)t 


Vqq (GeV) 


Cuu 


A 9 


canonical 


n fin 


U.Oo 








harder pair creation^ 1 ) 


0.60 


0.87 


~ 3 


~ 1 


~ 1/2 


softer pair creation^ 2 ) 


0.60 


0.43 


~ 1/2 


~ 1/3 


~ 1/6 


larger flux tube diameter® 


0.42 


0.58 


- 1/2 


- 1/2 


~ 1/2 



' recall /9 pc ~ 0.7 (3 pc so that the canonical value of (3 pc is ~ 0.4 GeV. 
« r g = 0.3 fm -> r g = 0.2 fm 
( 2 ) r q = 0.3 fm — > r g = 0.4 fm 
® / = 2 fm -> / = 1 fm 
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TABLE III. The fractions of the total nonresonant semileptonic decay rate predicted in low-lying 
exclusive channels. These fractions are obtained from Table I by taking (It^I 2 + 2|T 3/ / 2 | 2 )/A / 9 2 :ea 
(which may be obtained by summing the appropriate channels from Table I and dividing by 3) 
times flavor factors of | for / = 1, \ for I = t, and unity for K emission. Dominant subchannels 



can be read off from the 


single flavor rat< 


3S quote 


:d in Table I. 






percent of inclusive 


approximate phase space 


estimated net 


channel 


nonresonant rate in 


suppression factor in 


suppression after 




the heavy quark limit 


b — > c decays 


^■QCd/ttiq corrections 


(D + D*)7t 


0.9 




1/3 


0.5 


{D + D*) V + (£> + D*)rf 


0.3 




1/4 


0.4 


(D + D*)K 


0.6 




1/3 


0.5 


(D + D*)p 


2.8 




1/4 


0.4 


{D + D*)uj 


0.9 




1/4 


0.4 


(D + D*)K* 


1.8 




1/4 


0.4 


D T/2' K 


0.7 




1/4 


0.4 




0.2 




1/5 


0.3 




0.5 




1/4 


0.4 




0.7 




1/4 


0.4 


DT/2V + DT/2V' 


0.2 




1/5 


0.3 


DT1/2K 


0.5 




1/4 


0.4 


D T/2P 


2.2 




1/10 


0.15 




0.7 




1/10 


0.15 


DT3/2K* 


1.5 




1/10 


0.15 


DT/2P 


4.0 




1/5 


0.3 


D** /2 u 


1.3 




1/5 


0.3 


D*\ /2 K* 


2.7 




1/5 


0.3 


(D + D*)a 


1.7 




1/5 


0.3 
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(D + D*)f 


0.6 


1/5 


0.3 


(D + D*)K$ 


1.1 


1/5 


0.3 


(D + D*) ai 


2.8 


1/6 


0.25 


(D + D*)fi 


0.9 


1/6 


0.25 


(D + D*)K la 


1.9 


1/6 


0.25 


(D + D*)h 


2.4 


1/5 


0.3 


{D + D*)h! 


0.8 


1/5 


0.3 


{D + D*)K lb 


1.6 


1/5 


0.3 


{D + D*)a 2 


1.6 


1/6 


0.25 


(D + D*)f 2 


0.5 


1/6 


0.25 


(D + D*)K% 


1.1 


1/6 


0.25 



t For the separate emission of r\ and i]', or consideration of OZI-violating mixing in other 
nonets, mixing angles between the ui- and 0-like components may easily be introduced. Note 
that with the neglect of OZI mixing, emission of the 0-like meson is forbidden. 
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